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Abstract 

In this paper a computational method is introduced for proving that 
the diameter of the Rubik cube group must be less then or equal to 40 in 
the quarter turn metric. Recently it has been proven that when restricting 
the cube to edges the diameter is 18(see |3). Now given an arbitrarily 
scrambled cube one can solve the edges in at most 18 moves. When the 
edges are in place there are 44089920 possibilities for the corners. We 
prove that each of these states can be solved in at most 22 moves giving 
a new upper bound of 18-1-22=40. The program to solve this problem 
is made in the GAP language and takes about 10 hours to verify on a 
standard 3Ghz PC. Later some suggestions of further improvements on 
the diameter are given. 


General concepts 


In this section we introduce some general notions. The Rubik’s cube consists 
of 26 smaller cubies. We divide those cubies in three sets: the edge cubies ,the 
corner cubies and the center cubies. The edge cubies are those cubies that only 
have two colours and the corner cubies are those that have three colours. Finally 
the center cubies have only one color. So observe that an edge cubie has two 
faces and a corner cubie three faces. 


We will work with GAP. GAP is a computer algebra program and has a lot of 
special functions that makes handling of sets, permutations and other mathe¬ 
matical structures very easy. A manual with all instructions can be found at 
P]. When one wants to work with the cube group in GAP one has to specify 
its generators as permutations in cycle form. One suggestion of how this can be 
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done is given in Martin Schoenerts example Here the same example is used 
with a different enumeration: 
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32 16 39 

14 D 10 
41 7 34 


The picture above is a picture of the cube and we have enumerated each corner 
cubie face and each edge cubie face. U stands for up, L for left, R for right, D 
for down, F for front, B for back. 


We use the same letters for the generators as for the faces. For example the 
generator U twists the face ” up” clockwise. We define the generators by typing 
in GAP: 


U:=(l,2,4,8)(6,11,17,22)(25,26,29,36)(27,31,38,45)(33,40,46,44); 

L:=(3,6,12,20)(8,15,14,19)(25,27,32,42)(28,31,40,47)(35,41,46,36); 
F:=(4,9,16,12)(15,22,18,23)(27,33,43,35)(29,37,32,40)(36,45,39,47); 
R:=(2,5,10,18)(9,17,13,21)(26,30,39,33)(29,38,34,43)(37,45,44,48); 

B:=(1,3,7,13)(5,11,19,24)(25,28,34,44)(26,31,41,48)(30,38,46,42); 
D:=(7,14,16,10)(20,23,21,24)(28,35,37,30)(32,39,34,41)(42,47,43,48); 


Looking at these permutations one sees that the edges are numerated from 1 
to 24 and the corners from 25 to 48. We generate the cube group by typing in 
GAP: 


gap>cube:=Group(U,L,F,R,B,D); 

<permutation group with 6 generators> 

Now we want to let the cube act on the points 1 to 48. And we want the orbit 
of this action. This is done by typing in GAP: 


gap>Drbits(cube,[1..48],DnPoints); 

[[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24],[25, 
26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48]] 
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The response to the last command shows that points 1 to 24(edges) are permuted 
among themselves and the points 25 to 48(corners) are permuted among them¬ 
selves. This means that each group element of the cube can be uniquely written 
as the product of a permutation on the set 1 to 24 and a permutation on the set 
25 to 48. We will write it as ( 51 , 52 ) where gi is a permutation on the set 1 to 24 
and g 2 is a permutation on the set 25 to 48. The edge group is the group with 
generators U, L, F, B, R, D restricted to 1 to 24 i.e. if we write the generators 

as (Ug, Uc) — t/, /g) — T, {fej fc') — ipei ^c) — ^i ^c) — (^ej ^c) — ^ 

then the edge group is generated by Ue,le, fe,be,re, de- We denote it by Re- 
The corner group is defined to be the group with generators Uc, h, fc, be, Vc, dc- 
The Rubik cube group will be denoted by Gcube- Note that neither the edge 
group or the corner group is a subgroup of the cube. Actually if we for example 
intersect the corner group with the cube group we get a group two times smaller 
then corner group. The analogue holds for the edge group. This can be verified 
in GAP. The intersection of the corner group with the cube group we will call 
’’the group that hxes edges” or Ef and is very important in this text. 

The direct product between the corner group and the edge group gives a group 
twice the size of Rubik cube group and Rubik cube group is a subgroup of index 
2 of this direct product. This can be checked with GAP. 


Each element in a group can be written as an word where the generators of the 
group and the inverses of the generators are the letters. The length of an element 
of the group is defined as the minimal number of letters required to express it. 
The length of an element g is written L{g) when we use notation L{g) we mean 
the length of an element in Gcube with respect to generators {U, L, F, B, R, D}. 
And we write Le{g) when we mean the length of an element in the group Re 
with respect to generators {ue, h, fe, be, Ve, de}- 


The Symmetry Group 

Define M to be the group with generators 

fci = (1,24,16,22)(2,19,10,15)(3,21,12,17)(4,11,7,23)(5,14,18,8)(6,13,20,9) 

(25,48,32,45)(26,41,39,36)(27,44,42,37)(28,43,40,38)(29,31,34,47)(30,35,33,46) 
and 

k2 = (1,5,7,19)(2,10,14,8)(3,11,13,24)(4,18,16,15)(6,17,21,20)(9,23,12,22) 

(25,44,34,28)(26,48,41,31)(27,33,43,35)(29,37,32,40)(30,42,46,38)(36,45,39,47). 

This group acts on the set T = {U, L, F, B, R, D, U~^ , L~^, F~^, B~^, R~^, D~^} 
by conjugation i.e. if g G M and x € T then g~^xg € T. The action of conju¬ 
gation by ki can be written as a permutation with cycles: 

([/, B-\D, F-^){L, R-^){F, U-\B, D-^){R, R-^) 

and the action of ^ 2 : 
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(C/, i?, D, L)(C/-\i?-\£1-1, L-1). 

This action can be extended from T to the cube group. It is evident that two 
M conjugate elements have the same length. More details about this group can 
be found in pQ. 


Equivalence relations 

There are two important equivalence relations that we will refer to often in this 
text: 


1. Two elements gi, g 2 are equivalent if there 3h G M such that gi = h ^g 2 h. 
We refer to this relation as ~. 

2. Two elements 51,32 are equivalent if there 3h G M such that gi = h~^g 2 h 
or = h~^g 2 h . We refer to this relation as 

In both cases two equivalent elements have the same length. 
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Calculation of the diameter of the Edge Group 

(Re) 


We will also use that the diameter of the Cayley graph with respect to generators 
Ueje, fe,be, de^Te defined above is 18. This table is presented below and was 
found in The first column gives the number of elements of lengths 0 to 18. 
The second column gives the number of equivalence classes under relation ^ 
and the second column the number of equivalence classes under relation f«. 


Posted to Yahoo by Tom Rokicki on Jan 2, 2004 

Dist Positions Unique wrt M Unique wrt M+inv 


0 

1 

1 

1 

1 

12 

1 

1 

2 

114 

5 

5 

3 

1068 

25 

17 

4 

9819 

215 

128 

5 

89,392 

1,886 

986 

6 

807,000 

16,902 

8,652 

7 

7,209,384 

150,442 

75,740 

8 

63,624,107 

1,326,326 

665,398 

9 

552,158,812 

11,505,339 

5,759,523 

10 

4,643,963,023 

96,755,918 

48,408,203 

11 

36,003,343,336 

750,089,528 

375,164,394 

12 

208,075,583,101 

4,334,978,635 

2,167,999,621 

13 

441,790,281,226 

9,204,132,452 

4,603,365,303 

14 

277,713,627,518 

5,785,844,935 

2,894,003,596 

15 

12,144,555,140 

253,044,012 

126,739,897 

16 

23,716 

750 

677 

17 

30 

3 

3 

18 

1 

1 

1 


980,995,276,800 

20,437,847,376 

10,222,192,146 


This calculation implies that given an arbitrary element in the cube group we 
need to multiply it by at most 18 generators (in other words to apply at most 
18 moves to it) to take it in if/(the group that fixes edges). 
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Mapping the corner group on numbers 


It is well known that the corner group is isomorphic to S'g k Zg. This means 
physically that each configuration in the corner group can be described as a 
permutation of the cubies and an orientation of the cubies. Let (j) : i—>- 

Ss K Zg be such a isomorphism. If g G Rc then (j){g) can be uniquely written 
as 4 'i 9 )i 4 '{ 9 ) 2 ], (pig)! G S 8 , 4 '{g )2 s it is obvious how to find a bijection 
from Zg to [0,3"^ — 1]. Call this bijection ipi. There is a simple way to define 
a bijective map 'il’2 '■ Ss [0,8! — 1]. One can define a natural order relation 
on permutations in Ss which means that given two permutations ai , (J2 and 
(Ti < (72 then either cti(1) < (72(1) or if n is the biggest positive integer with the 
property ai{k) = a 2 {k),k = 1, ..,n then cti( n +1) < a 2 {n + l). We simply define 
V' 2 (ct) = S{t|t < cr}. Now simply define a bijection f : R^ ^ [1,88179840] as 
f ig) = • ip2{(j){g)i) + ■fi{(j){g)2)- 


The main objective of this paper 


In this paper we show how to construct explicit expressions in generators of 
length 22 or less for each element in if/(’’the group that fixes edges”). This 
group has 44089920 elements, this can be shown by typing in GAP: 


gap> cube:=Group(U,L,F,B,R,D); 
<permutation group with 6 generators> 
gap> corners:=Group(uc,lc,fc,bc,rc,dc); 
<permutation group with 6 generators> 
gap> Size(Intersection(corners,cube)); 
44089920 


Since we need at most 18 moves to take an element to Ej and at most 22 moves 
to take an element of A/ to the identity. This means also that the diameter of 
the Rubik’s cube group is less then or equal to 40. 


Algorithm to express elements of Ef in generators 


We need to introduce some notations. Define the following function p : Gcube 
Re by p((( 7 i, 32 )) = 5 i-We call this map ’’projection on the first component”. 
Now for a fix g G define the set S(n,g) = {ig,x) G Gcube\L{{g,x)) < n} 
i.e. the set of all g in the cube group that can be expressed as a product 
of at most n generators and p(g) = g for all g G S{n,g). One can easily 
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see that S{n,g) is empty if Le{g) > n. Define S(k,g) ■ S(n,h) = {x ■ y\x € 
S{k,g),y G S{n, h)}. It is obvious that >S'(10, g) ■ <S'(12, g~^) = S{22, id). 

So Ug6ii^,L^(g)<io ‘S'(10, g) ■ >S'(12, g ^) = S(22, id). 

Following relations are helpful when trying to understand the algorithm we will 
describe: 


1. h ^S{g,n)h = S{h ^gh,n),h€M 

2. h~^S{k, g) ■ S{n, g~^)h = S{k, h~^gh) ■ S{n,h~^g~^h),h G M. 

3. {S{g,n))-^ = S{g-'^,n) 


In order to keep track on which elements we have found generator expressions 
for, a ’’bit vector” is used. A ’’bit vector” is a vector in which each entry can take 
only two values 0 or 1. GAP uses ’’false” for 0 and ’’true” for 1. Initially all en¬ 
tries in this ’’bit vector” are set to 1 or ’’true”. A bijection f : Rc [1, 88179840] 
was defined on page 6. When generator expressions for a given element g G Ej 
have been found entry number f{g) is set to 0 or ’’false”. When we use the 
expression ’’element g has been checked of’ we mean that entry f{g) in the ’’bit 
vector” has been set to ” false”. 

We also need to collect information about how a given element can be expressed 
in generators. In GAP we use for this a set we call ’’positions”. The algorithm 
we describe stores only one representative from each equivalence class under 
relation «. When we write ’’store generator expression” we mean that the gen¬ 
erator expression will be saved to the set ” positions”. 


Now we are ready to describe the main step of the algorithm: 

Define A to be the set A = {g Re\Le{g) < 10}. Fix a g G A and do following 
steps: 


1. Compute S{10,g) and generator expressions for its elements 

2. Compute S'(12,g“^) and generator expressions for its elements 

3. Compute S{10,g) ■ S{12,g~^) and generator expressions for its elements 

4. For each element elm in ^(lO,^) • S{12,g~^) do following steps: 

(a) Check if elm has been checked. If yes do nothing if no do the following 
steps 

i. Store generator expression of elm 

ii. Compute the orbit of elm under relation «. 

iii. Check of every element in the orbit. 

5. Compute S'(10,5“^) • S{12,g) and generator expressions for its elements 
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6 . For each element elm in 5'(10, g • S{12,g) do following steps: 

(a) Check if elm has been checked. If yes do nothing if no do the following 
steps 

i. Store generator expression of elm. 

ii. Compute the orbit of elm under relation w. 

iii. Check of every element in the orbit. 

Let 0{g) be the orbit of g under relation w. Then the above algorithm checks 
of every element in £'( 10 , 5 ) • S{12,g-^). 

The algorithm consists of repeating above procedure for enough elements g € A 
which are not equivalent under relation «. It should be noted that building 
£(12, g) is the most time consuming step in the algorithm. Once we have £(12, g) 
it takes no time to build £(12, h~^gh). Thats why we do the above steps for 
the whole equivalence class. 

If we for example run the algorithm with g choused as identity than 3079007 
elements are checked of. 


Optimal solvers 


As mentioned above £( 22 , id) can be written using following formula: 

U £(10, g). £(12,5-1)= £(22, zd) 

geRe,L.,{g)<10 

In the above algorithm we see how to use it in order to compute each element 
in Ef as a product of generators. Obviously some elements we compute might 
have length that is less then 22. To compute those we use the same algorithm 
with the exception that £( 10 , 5 ) is replaced by £( 8 , 5 ). And £(20, zd) can be 
written as: 

IJ £( 8 , 5 )-£( 12 , 5 "') = £( 20 ,zd) 

gGRe ,-^e (5) 

Now we have all elements of length less then or equal to 20 in Ay. To get only 
elements of length exactly 22 we must take the difference between £( 22 , id) and 
£(20, zd). We continue computing the set £(18, zd) using the formula: 

U £( 8 , 5 )-£( 10 , 5 -') = £(18,zd) 

geRe,Le(g)<8 


And finally we compute £(16, id), £(14, zd), £(12, zd), £(10, id) and £( 8 , id). To 
compute elements of length n we simply take the difference £(n, 5 ) — £(n — 1 , 5 ). 
The program for this calculation is not given here but we give the result on next 
page. Observe that 22q means length 22. 



Distance 

Nr of pos 

Unique wrt 

M Unique wrt M + inv 

Oq 

1 

1 

1 

2q 

0 

0 

0 

4q 

0 

0 

0 

6q 

0 

0 

0 

8q 

240 

5 

3 

lOq 

288 

6 

3 

12q 

8764 

197 

113 

14q 

116608 

2475 

1303 

16q 

840513 

17620 

9091 

18q 

9884342 

206424 

104523 

20q 

30623418 

639362 

323859 

22q 

2615746 

54895 

28528 


44089920 

920980 

467424 

We have also computed the length of all elements in the ’’Group that fixes 

cubies”. This 

is the set of all elements in the Rubik cube group such that the 

cubies are in 

their correct spot but their orientation is arbitrary. This is an 

abelian normal subgroup of the cube. This 

can be checked with GAP. The 

order of this group is 4478976. 

Here is the table with lengths of each element: 

Distance 

Nr of pos 

Unique wrt 

M Unique wrt M + inv 

Oq 

1 

1 

1 

2q 

0 

0 

0 

4q 

0 

0 

0 

6q 

0 

0 

0 

8q 

0 

0 

0 

lOq 

0 

0 

0 

12q 

441 

11 

8 

14q 

3944 

87 

52 

16q 

32110 

708 

396 

18q 

456025 

9656 

5009 

20q 

2873194 

60399 

30978 

22q 

1113236 

23652 

12424 

24q 

25 

4 

4 


4478976 

94518 

48872 


One of the positions of length 24 is the well known ’’superflip”. This position 
alone generates the center of the Rubik cube group. The other positions up to 
conjugation with M are: 
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F R’ U B2 U R D F U’ B R B2 D B R’ U’ F' B R U2 D’ (24q*,21f) 

F R’ F2 R’ U F U L B2 R B D2 F R' B’ D B L’ D’ R’ U’ (24q*, 21f) 

F R D2 B’ R’ F R’ D2 B' D’ L F L2 U B’ U F L2 U' B’ (24q*, 20f) 


This positions were the only ones not computable in 22q and the expressions in 
generators were found using Michael Reids solver see [21- 


Upperbounds on different subsets 


Take an element g G Re and L{g) = n {L{g) is the length of g) in Re. Then for 
any element {g,y) G Gcube we need at most n moves to make it an element of 
the form {id, x) G Gcube-And an element {id, x) G Gcube needs at most 22 moves 
applied on it to make it the identity. So L{{g, y)) < 22 + n in Gcube-There is only 
one element h in Re with L{h) = 18 in i?e(according to the above diagram). 
So the only elements in Gcube that may need 40 moves are those in the set 
S = {{h, x) G Gcube}-If a solver solves each position in S in less then 40 moves 
then one gets a new upper bound 39. Actually one only needs to solve a position 
from each equivalence class under relation « which gives about 500000 positions 
needed to be solved. This computation has not been done here. But soon after 
the author posted generator expressions for elements in i?/ this has been done 
by another author. See [HI for details. 


The program 


To begin with we are going to create some sets we are going to work with. The 
sets will be presented in GAP as lists. The first one is s. 

This set is made by typing in gap: 


gap>s;=[U,L,F,B,R,D,U'-l,L'-l,F~-l,B~-l,R~-l,D~-l]; 


The other sets are: 


Some sets 

veven the set of all elements in the edge group of length 0, 2, 4, 6. 
vodd the set of all elements in the edge group of length 1, 3, 5. 
vevenr the set of all elements in the Rubik cube group of length 0, 2, 4, 6. 
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voddr the set of all elements in the Rubik cube group of length 1, 3, 5. 

The elements of this sets are stored as permutations with disjoint cycles. 

vevenc Contains at entry k the same element that vevenr contains at entry 
k but stored as products of generators. Throughout this program we store a 
product of generators as a vector whose entries are integers between 1 and 12. 
Integer n means generator s[n] and s was defined above, 
voddc is analogue to vevenc but the set vevenr is replaced with voddr. 

All the sets containing permutations with cycles are sorted with respect to the 
natural order relation on permutations which means that given two permuta¬ 
tions CTi, (72 and (Ti < (72 then either cti(1 ) < (72(1) or if n is the biggest positive 
integer with the property cri(fc) = 02 {k), k = 1 ,.., n then ai{n + 1) < a 2 {n + 1). 
The elements in the set vevenr are ordered in following way: 
{gi,hkj,{gi,hk 2 ),-,igi,hk^),{g 2 ,hi^),(g 2 ,hi^),.... This is because the first 
component is a permutation on the set 1 to 24 while the second component 
a permutation on the set 25 to 48. Then the set veven is ordered in following 
way: 

gi, g 2 , ■■■ Note that we refer to an element of vevenr as an element of the form 
{g, h) although it is stored as a permutation with disjoint cycles. 

We use two more vectors yl and y2 where the first has the same size as vodd 
and the second the same size as veven.The vector yl contains at position k a 
number that gives the position in voddr where we have the first occurrence of 
an element of the form (vodd(fc), h). 

Finally we call the ’’bit vector” for rev which has the same length as Rc 
(88179840 elements). When the program starts all its entries are set to true. 


Some functions 

We will predefine some functions in order to make the program code easier: 


MakeSetv(n) if n=l then it returns set [vevenr, voddr] if n=2 returns set 
[veven, voddj. 

MakeSetvc() returns [vevenc,voddc]. 

MakeSety(n) if n=l returns set yl if n=2 returns set y2. 

MakeSet(n, g) This function will give a set of the form {g,h),h € Rc and all its 
elements of length n or less. The returned set is of the form [setl,set2] where 
the elements of setl are expressed as the elements of vevenr and the elements 
of set2 are expressed as products of generators. 

PerniToNr(g) This is a bijective function PermToNr : Rc 1 -^ [1,88179840]. 
EqClassRel2(el) Returns the set {PermToNr{g)\g « el}. 

WriteRcv(setl) For each x S setl does rcv{x) = false. 

Redf(setl).Partitions set setl into equivalence classes under relation w defined 
above and returns a set containing a representative from each class. 
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CharToPerm(x) is a function that converts a text string to a permutation. 
This was done so that the attached data don’t take to much space in this paper. 
MultEl(x) if for example x := [1,2,1,4,11, 2] then this function returns s[l] * 
s[2] * s[l] > 1 = s[4] * s[ll] s[2]. 

Invg(x) if for example x := [1,2,1,4,11,2] then it returns [8,5,10,7,8,7] i.e. 
MultEl(x)-i = MultEl(Invg(x)). 
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The Code 


Before we describe the code some GAP commandos must be explained. Those 
can also be found in the manual |3] but we describe them here for completeness. 
To build a set in GAP one types: 


gap>set: = [] ; 


To build a set of sets type: 


gap>set : = [[],[],[],[]]; 


The above set contains 4 empty subsets. To store the identity in the hrst subset 
type in GAP: 


gap>set[1] [1]: = (); 

Now the first set in the set ’’set” contains the identity. 

To build a ’’bit vector” of size 1000 with all entries set to ’’true” or 1 just type: 
gap>vec:=BlistList([1..1000],[1..1000]); 

To build a ’’bit vector” with all entries false just type: 
gap>vec:=BlistList([1..1000],[]); 

If we type in GAP: 


gap>set:=[1,1]; 


then we have build a set that contains the same element on two different entries, 
maybe it is better to call it a list. But if we type: 


gap>set2:=Union(set,set); 

then the set ”set2” will contain only element 1. If we instead of union would have 
taken intersection the same thing should have happened. A union, intersection 
or difference of two sets will never return a set with doubles. Furthermore after 
performing those operations the set will even be sorted. 
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Generators and the Symmetry Group 


We define all generators we are going to use and the symmetry group in following 
file: 


#Generators for the cube group 

U:=(l,2,4,8)(6,11,17,22)(25,26,29,36)(27,31,38,45)(33,40,46,44); 

L:=(3,6,12,20)(8,15,14,19)(25,27,32,42)(28,31,40,47)(35,41,46,36); 

F:=(4,9,16,12)(15,22,18,23)(27,33,43,35)(29,37,32,40)(36,45,39,47); 
R:=(2,5,10,18)(9,17,13,21)(26,30,39,33)(29,38,34,43)(37,45,44,48); 

B:=(1,3,7,13)(5,11,19,24)(25,28,34,44)(26,31,41,48)(30,38,46,42); 
D:=(7,14,16,10)(20,23,21,24)(28,35,37,30)(32,39,34,41)(42,47,43,48); 

#Generators for edge group an corner group 

ue:=RestrictedPerm(U,[1..24]);uc:=RestrictedPerm(U,[25..48]); 
le:=RestrictedPerm(L,[1..24]);Ic:=RestrictedPerm(L,[25..48]); 
fe:=RestrictedPerm(F,[1..24]);fc:=RestrictedPerm(F,[25..48]); 
re:=RestrictedPerm(R,[1..24]);rc:=RestrictedPerm(R,[25..48]); 
be:=RestrictedPerm(B,[1..24]);bc:=RestrictedPerm(B,[25..48]); 
de:=RestrictedPerm(D,[1..24]);dc:=RestrictedPerm(D,[25..48]); 

#Generators for the group of symmetries 

kl:=(1,24,16,22)(2,19,10,15)(3,21,12,17)(4,11,7,23)(5,14,18,8)(6,13,20,9) 
(25,48,32,45)(26,41,39,36)(27,44,42,37)(28,43,40,38)(29,31,34,47)(30,35,33,46); 
k2:=(1,5,7,19)(2,10,14,8)(3,11,13,24)(4,18,16,15)(6,17,21,20)(9,23,12,22) 
(25,44,34,28)(26,48,41,31)(27,33,43,35)(29,37,32,40)(30,42,46,38)(36,45,39,47); 
#Generators for the group of symmetries restricted to edge resp corners, 
kle:=RestrictedPerm(kl,[1..24]);k2e:=RestrictedPerm(k2,[1..24]); 
klc:=RestrictedPerm(kl,[25..48]);k2c:=RestrictedPerm(k2,[25..48]); 

#symmetry group 
M:=Group(kl,k2); 

#symmetry group restricted to corners 
Mc:=Group(klc,k2c); 

#symmetry group restricted to edges 
Me:=Group(kle,k2e); 


This code will be saved to file gens.txt. 
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Function MakeSetv 


This function builds the sets veven, vodd, vevenr, voddr. The function hrst 
selects generators depending on if we want to build veven, vodd or vevenr, 
voddr. If we want to build vevenr, voddr we type: 


gap>a:=MakeSetv(1); 


The function will return a set of sets containing the sets vevenr and voddr. In 
the above case a[l] is vevenr and a[2] is voddr. Lets describe the algorithm that 
this program uses: 

1. First a set of sets is build containing 6 empty sets and a set containing 
the identity: 

v: = [[()], [],[],[],[], [],[],[]] 

2. Now take each element v[l] and multiply it be each generator and store 
the result in the next set i.e. v[2]. 

3. Now take v[2] and multiply each element of it with the generators and 
store the result in v[3]. Remove doubles from v[3] and elements contained 
in v[l]. 

4. Now take the next set in v i.e. v[3] multiply each element by generators, 
store the result in v[4] remove doubles from v[4] and elements contained 
in v[2]. 

5. Continue doing this procedure until v[6] i.e. v[6] is the last set of which 
elements will me multiplied by generators and the result stored in v[7] 

6 . Now return the set {Union(v[l],v[3],v[5],v[7]),Union(v[2],v[4],v[6])} i.e. {veven,vodd}. 
See next page for the code of above described function: 
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MakeSetv:=function(n) 

local veven,vodd,vevenr,voddr,gens,selct,v,k,x,z,y,cnt; 
gens: = [[()], [()]];selct:=n; 

gens[l]:=[U,L,F,B,R,D,U~-1,L~-1,F~-1,B~-1,R~-1,D~-1]; 
gens[2]: = [ue,le,fe,be,re,de,ue"-l,le"-l,fe~-l,be~-l,re"-l,de"-l] ; 
v: = [[()],[()],[()],[()],[()],[()],[()]]; 
f or y in [1. . 6] do 
k: = 

cnt:=1; 

for X in v[y] do 

for z in [1..12] do 

k[z] [cnt]:=x*gens [selct] [z]; 

od; 

cnt:=cnt+l; 

od; 

v[y+l] :=Union(k[l] ,k[2] ,k[3] ,k[4] ,k[5] ,k[6] ,k[7] ,k[8] ,k[9] ,k[10] ,k[ll] ,k[12]) ; 
if y>l then v[y+1];=Difference(v[y+1],v[y-1]);fi; 

od; 

if selct=l then 

vevenr:=Union(v[1],v[3],v[5],v[7]);voddr:=Union(v[2],v[4],v[6] ); 
return [vevenr,voddr]; 
elif selct=2 then 

veven :=Union(v[l] ,v[3] ,v[5] ,v[7] ) ; vodd: =Union(v [2] ,v[4] ,v[6] ) ; 
return [veven,vodd]; 

fi; 

end; ; 

This code will be saved to file makesetv.txt. 
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Function MakeSetvc 


This function will build the sets vevenc and voddc. We remind that this sets 
contain the same number of elements as vevenr and voddr. We have defined 
a list ”s” above as s = [U,L,F,B,R,D,U-^,L-^,F-^,B -^To 
define it in GAP just type: 

gap>s: = [U,L,F,B,R,D,U'-l,L"-l,F"-l,B"-l,R"-l,D"-l] ; 

We now want to store generator expressions for elements in vevenr and voddr 
in sets vevenc and voddc if for instance vevenc contains at position 10 generator 
expression [1,2,7] then vevenr should contain at position 10 element s[l]*s[2]*s[7]. 
We will describe an algorithm how to achieve this. But hrst we need a few 
preparations. GAP has a function called ’’Tuples” that we will use. We explain 
it by giving an example: 


gap>Tuples([1..5],2); 


After typing this GAP will return a set containing all possible elements of form 
(x, y) with X G {1, 2,3,4, 5}, y G {1, 2,3,4, 5}. And we also use a function called 
PositionSet. We explain it using an example: 


gap>PositionSet(veven,(1,2,3)); 


If (1,2,3) is contained in veven then this function will return a number n such 
that veven[n] = (l,2,3). If it is not contained it will just return fail. Now we are 
ready to describe the algorithm. It is described only for vevenc since voddc is 
analogue: 

1. Build a ’’bit vector” with same size a vevenr with all elements set to 
’’false”. 

2. Build the set of all tuples on elements {1,2,3,4,5,6,7,8,9,10,11,12} and 
length 2. Call this set Tu2. 

3. Take each element {x,y) G Tu2 and Hnd the position in vevenr of the 
element s[x] * s[y]. Call this position n. If entry number n is ’’false” in the 
’’bit vector”, set it to ’’true” and put element {x,y) in vevenc at position 
n. If entry number n in the ’’bit vector” is ’’true” then do nothing. 

4. Repeat step 2 and 3 for tuples of length 4 and 6 on the set {1,2,3,4,5,6,7,8,9,10,11,12}. 
The code for this function is given on next page. 
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s: = [U,L,F,B,R,D,U'-l,L"-l,F"-l,B~-l,R~-l,D~-l] ; 

MakeSetvc:=function() 

local binv,binv2,vevenc,voddc,x,set,z,elm,g,pos; 

binv:=BlistList([1..Size(vevenr)],[]); 

binv2:=BlistList([1..Size(voddr)],[]); 

vevenc : = [] ; 

voddc : = [] ; 

for X in [1..6] do 

if x=l or x=3 or x=5 then 

set:=Tuples( [1..12],x); 

for z in set do 

elm:=(); 

for g in [l..Size(z)] do 

elm:=elm*s[z[g]]; 

od; 

pos:=PositionSet(voddr,elm); 

if binv2[pos]=false then voddc[pos]:=z;binv2[pos]:=true;fi 
od; 

fi; 


if x=2 or x=4 or x=6 then 
set:=Tuples( [1..12],x); 
for z in set do 
elm:=(); 

for g in [l..Size(z)] do 

elm:=elm*s[z[g]]; 

od; 

pos:=PositionSet(vevenr,elm); 

if binv[pos]=false then vevenc[pos]:=z;binv [pos]:=true;fi; 
od; 

fi; 

od; 

vevenc [1] : = [] ; 
return[vevenc,voddc]; 
end; ; 
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Function MakeSety 


As described before vector yl has the same size as vodd and vector y2 the same 
size as veven. We only describe how to build yl since y2 is analogue. Remember 
that if an element g G Re has position n in veven then y[n] contains a number 
m that gives the position in vevenr of the first element who’s projection on Re 
is g. By projection we mean the map (g, x) g. We use following algorithm to 
build y2: 

1. Build a vector ’’vtemp” the same size as vevenr that contains at position 
n the projection of vevenr[n] on Re- 

2. Do following for each g G veven. 

(a) Find position of g in veven call it ”pos” 

(b) Find the position of the first occurrence of g in vtemp call it ”pos2”. 

(c) Set y2[pos]=pos2 

The code is described below: 


MakeSety:=function(n) local vtemp,v,vr,x,y,cnt,g; 

vtemp:=[()];cnt:=l; 

if n=l then v:=vodd;vr:=voddr; 

elif n=2 then v:=veven;vr:=vevenr; 

fi; 

x:=[l..Size(vr)];y:=[l..Size(v)]; 
for g in X do 

vtemp[g]:=RestrictedPerm(vr[g],[1..24]); 

od; 

for g in V do 

y[PositionSorted(v,g)]:=PositionSorted(vtemp,g); 

od; 

return y; end;; 


This code will be saved to file makesety.txt 


Function PermToNr 

Here is the program described in the section ”How to map a permutation on a 
number”. It takes a g G Re and maps it on S'g x Zg. Lets describe how to map 
a permutation of Ss on numbers [1,8!] in more detail. Given a a G Ss we want 
to determine all elements that are less then a with respect to the order relation 
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that was defined on page 11. This is a simple combinatorial problem and we 
give an example of how to determine the number of elements less then a given 
permutation and then generalize this to an arbitrary one. Take for example 


1 2 3 4 5 6 7 8 

5 6 3 4 2 1 8 7 


We divide the set P of permutations less 


than cr in following disjoint subsets: 


1. The subset with {t(1) < 5,r G P}. Then t(1) can be any element in the 
set {1,2,3,4} so this set has size 7! • 4 

2. The subset with {t(1) = 5,r G P}. Then t(2) can be any element in the 
set {1,2,3,4} so the size of this set is 6! • 4. 

3. The subset with {t(1) = 5, r(2) = 6, t G P} then r(3) can be any element 
in the set {1, 2} and the size of this set is 5! • 2. 

4. The subset with {t(1) = 5,t(2) = 6,t(3) = 3,t G P} then r(4) can be 
any element in the set {2,1}. The size of this subset is 4! • 2. 

5. The subset with {t(1) = 5,t(2) = 6, r(3) = 3, r(4) = 4} then t(5),t G P 
is in the set {1}. The size of this set is 3!. 

6 . The subset with {t(1) = 5, r(2) = 6, r(3) = 3,t(4) = 4, t(5) = 2,r G P} 
and t(6) can be any element in the set {}. The size of this subset is 0. 

7. The subset with {r(l) = 5, r(2) = 6,t(3) = 3,r(4) = 4,r(5) = 2,t(6) = 
l,r G P} then t(7) can be any element in the set {7}. The size of this 
subset is 1!. 


Since all of this subsets are disjoint and their union is P it follows that the size 
ofPis4-7! + 4-6! + 2- 5! + 2-4! + 3! + 0 + l!. This can be easily generalized for an 

1 2 3 4 5 6 7 8 

cr(l) cr(2) ct(3) ct(4) (t( 5) cr(6) cr(7) cr(8) 

then the number of permutations less than a is 


arbitrary cr. Let cr = 


^(7 - fc + 1)! • {a{k) - 1) • tt{CT(l) < aik) | / G {1,.., k - 1}} 

fc=i 

The code for bijection f ■. Rc^ [1,88179840] can be seen on next page. 
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PermToNr:=function(g) local ml,m2,m3,m4,m5,m6,m7,m8,y,w,nr,kl,k2; 
y:=[l..48]; 

y[25] :=l;y[31]:=2;y[46]:=3;y[26]:=l;y[38]:=2;y[44]:=3;y[27]:=l;y[36]:=2;y[40]:=3; 
y[28]:=l;y[41]:=2;y[42]:=3;y[29]:=l;y[33]:=2;y[45]:=3;y[30]:=l;y[34]:=2;y[48]:=3; 
y[32] :=l;y[35]:=2;y[47]:=3;y[37]:=l;y[39]:=2;y[43]:=3; w: = [l..48]; 
w[25]:=l;w[31]:=l;w[46]:=l;w[26]:=2;w[38]:=2;w[44]:=2;w[27]:=3;w[36]:=3;w[40]:=3; 
w[28]:=4;w[41]:=4;w[42]:=4;w[29]:=5;w[33]:=5;w[45]:=5;w[30]:=6;w[34]:=6;w[48] :=6; 
w[32] :=7;w[35] :=7;w[47] :=7;w[37] :=8;w[39] :=8;w[43] :=8; 

ml:=w[25~g];m2:=w[26"g];m3:=w[27~g];m4:=w[28~g];m5:=w[29"g];m6:=w[30~g];m7:=w[32~g]; 

m8:=w[37~g]; if ml<m2 then m2:=m2-l; fi; if ml<m3 then m3:=m3-l; fi; 

if m2<m3 then m3:=m3-l; fi; if ml<m4 then m4:=m4-l; fi; if m2<m4 

then m4:=m4-l; fi; if m3<m4 then m4:=m4-l; fi; if ml<m5 then 

m5:=m5-l; fi; if m2<m5 then m5:=m5-l; fi; if m3<m5 then m5:=m5-l; 

fi; if m4<m5 then m5:=m5-l; fi; if ml<m6 then m6:=m6-l; fi; if m2<m6 

then m6:=m6-l; fi; if m3<m6 then m6:=m6-l; fi; if m4<m6 then 

m6:=m6-l; fi; if m5<m6 then m6:=m6-l; fi; if ml<m7 then m7:=m7-l; 

fi; if m2<m7 then m7:=m7-l; fi; if m3<m7 then m7:=m7-l; fi; if m4<m7 

then m7:=m7-l; fi; if m5<m7 then m7:=m7-l; fi; if m6<m7 then 

m7:=m7-l; fi; 

kl:=(ml-l)*7*6*5*4*3*2+(m2-l)*6*5*4*3*2+(m3-l)*5*4*3*2+(m4-l)*4*3*2+ 

(m5-l)*3*2+(m6-l)*2+m7-l; 

k2:=(y[25-g]-l)*3-0+(y[26-g]-l)*3~l+(y[27-g]-l)*3-2+(y[28-g]-l)*3-3+ 
(y[29-g]-l)*3~4+(y[30 "g]-1)*3"5+(y[32~g]-1)*3"6; nr:=2187*kl+k2+l; 
return nr; end;; 


This code will be saved to file permtonr.txt. 

Nobody is expected to follow this code. That is why a program is given to verify 
that this function is a bijection and it is all that is needed for the proof: 


corners:=Group(uc,Ic,fc,bc,rc,dc); rev:=BlistList( [1..88179840], []); 
for X in corners do 
rev[PermToNr(x)]:=true; 
od; 

Print(rcv=BlistList([l..88179840],[1..88179840]),"\n"); 


If the program returns true then it is a surjection and consequently a bijection. 
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The function MakeSet(n,g) 

This function builds S'(10,g) and S'(12,g). To build 5'(10,(/) just type: 
gap>MakeSet(10,g); 

We explain the algorithm for S'(12,g) 

1. Build an empty set S 

2. For each element h in veven do 

(a) Build the set Sh = {{h, x) S vevenr} and the set Sh-^g = {{h~^g, x) G 
vevenr}. 

(b) Build the set {gi • 32 | G Sh,g 2 G S'/t-ig} and put each element 
from this set in S. 

3. Return the set S (or the set S = S{12,g) 

The implementation of this function as code can be seen on next page. 
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MakeSet:=function(n,g) local set,y,veo,veor,a,b,z,x,kl,j,k,l,m,t,set2,set3,nre,veoc 
if n=10 then 

y;=y1;veo;=vodd;veor:=voddr;veoc:=voddc; 
elif n=12 then 

y;=y2;veo:=veven;veor:=vevenr;veoc:=vevenc; 

fi; 

set:=[];t:=[1..Size(veo)];kl:=Size(veo);set2:=[];set3:=[];nre:=0; 
for X in t do 
if x=kl then 

a:=Size(veor)+l-y[x]; 
else 

a:=y[x+l]-y[x] ; 

fi; 

z:=PositionSet(veo,veo[x]~-l*g); 
if zOfail then 

if z=Size(veo) then b:=Size(veor)+l-y[z]; 
else b:=y[z+1]-y[z]; 

fi; 

1: = [1. .a] ;m: = [l. .b] ; 


for j in m do 
for k in 1 do 

if PositionSet(set,veor[y[x]+k-l]*veor[y[z]+j-l])=fail then nre:=nre+l; 

AddSet(set,veor[y[x]+k-l]*veor [y[z]+j-l]); 

set2[nre]:=veor[y[x]+k-l]*veor[y[z]+j-l] ; 

set3[nre]:=Concatenation(veoc[y[x]+k-l],veoc[y[z]+j-l]); 

fi; 

od; 

od; 

fi; 

od; 

return [set2,set3];end;; 


This code will be saved to file makeset.txt 


23 



The function EqClassRel2(el) 


This function returns a set that contains the image under the bijection f : Rc 
[1, 88179840] of every element equivalent to el under relation «. The algorithm 
is evident and need not to be explained. 


EqClassRel2:=function(g) local h,set; 

set: = [] ; 

for h in Me do 

AddSet(set,PermToNr(h~-l*g*h)); 
AddSet(set,PermToNr(h~-l*g~-l*h)); 

od; 

return set; end;; 


This code will be saved to file eqclassrel2.txt 


The function WriteRcv(set) 


This function has as input a set of indices or numbers and is called by command: 


gap>WriteRcv(set); 


The command above takes every n G set and sets entry n in the ’’bit vector” to 
’’false” or 0. 


WriteRcv:=function(set) 
local x; 
for X in set do 
rev[x]:=false; 
od; 
end; ; 

This code will be saved to file writercv.txt 
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The function Redf(set) 


As mentioned earlier in the text the algorithm finds generator expressions for 
elements in A/ and each time it does this we need to give it a g G Re- It was 
also mentioned that if given the algorithm a, g G Re that is equivalent under 
relation « with another h € Re that has been before, then this h & Re will not 
yield any new generator expressions. Redf{set) returns a set containing only 
one representative from each equivalence class under relation w. The algorithm 
is simple and can be seen directly by looking at the code. 


Redf:=function(set) local x,y,set2,setemp,vec; 
vec:=BlistList( [1..Size(set)],[1..Size(set)]);set2: = [] ; 
for X in [1..Size(set)] do 
setemp : = [] ; 
if vec[x] then 
for y in Me do 

vec[PositionSet(set,y~-l*set[x]*y)]:=false; 
vec[PositionSet(set,y~-l*set [x]~-l*y)]:=false; 

AddSet(setemp,y~-l*set[x]*y); 

AddSet(setemp,y~-l*set [x]~-l*y); 
od; 

if lsEmpty(lntersection(setemp,set2)) then AddSet(set2,set[x]); fi 

fi; 

od; 

return set2; end;; 


This code will be saved to file redf.txt 
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The function CharToPerm(x) 


Part of the g € Re we are giving as input to the algorithm will come from 
the set Redf {veven). With this set we find generator expression for 43984820 
elements in Ef. To find expressions for the rest of 105100 we need some more 
g G Re that has been extracted from a large list and are given at the end of 
this paper. This list took about 2 days to find on 20 computers. It consisted 
of representatives under relation Ri of elements of the set {g G Re,L{g) = 8} 
and very few from {g G Re, L{g) = 10}. From the last set only 17. Those were 
found using a solver found at |2]. In order to compress the space the elements 
have been well coded and in order to make them permutations in Re we need 
the function CharToPerm{x). If we want to convert an element of this list to 
a permutation we type in GAP: 


gap>CharToPerm("acekrihownvs"); 
(4,6,5,21,23,15,7,22,8,13,10,16,12,24)(14,20) 


Since this is just something technical it needs no further explanation. 

CharToPerm:=function(g) 
local s,v,w,m,k; 
v:=[l..24]; 

v[l]:=ll;v[2]:=17;v[3]:=19;v[4]:=22;v[5]:=13;v[6]:=8;v[7]:=24;v[8]:=6;v[9]:=18; 
v[10]:=21;v[ll]:=l;v[12]:=15;v[13]:=5;v[14]:=20;v[15]:=12;v[16]:=23;v[17]:=2; 
v[18]:=9;v[19]:=3;v[20]:=14;v[21]:=10;v[22]:=4;v[23]:=16;v[24]:=7; 
w:=[l..256];k:=lnputTextString("abcdefghijklmnopqrstuvwx"); 
w[ReadByte(k)]:=1;w[ReadByte(k)]:=ll;w[ReadByte(k)]:=2;w[ReadByte(k)]:=17; 
w[ReadByte(k)]:=3;w[ReadByte(k)]:=19;w[ReadByte(k)]:=4;w[ReadByte(k)]:=22; 
w[ReadByte(k)]:=5;w[ReadByte(k)]:=13;w[ReadByte(k)]:=6;w[ReadByte(k)] :=8; 
w[ReadByte(k)]:=7;w[ReadByte(k)]:=24;w[ReadByte(k)]:=9;w[ReadByte(k)]:=18; 
w[ReadByte(k)]:=10;w[ReadByte(k)]:=21;w[ReadByte(k)]:=12;w[ReadByte(k)] :=15; 
w[ReadByte(k)]:=14;w[ReadByte(k)]:=20;w[ReadByte(k)]:=16;w[ReadByte(k)]:=23; 
m:=lnputTextString(g);s:=[1..24]; 

s[l] :=w[ReadByte(m)] ;s[v[l]] :=v[s[l]] ; s [2] :=w[ReadByte(m)] ;s[v[2]] :=v[s[2]] ; 
s [3] : =w [ReadByte (m) ] ; s [v [3] ] : =v [s [3] ] ; s [4] : =w [ReadByte (m) ] ; s [v [4] ] : =v [s [4] ] ; 
s [5] : =w [ReadByte (m) ] ; s [v [5] ] : =v [s [5] ] ; s [6] : =w [ReadByte (m) ] ; s [v [6] ] : =v [s [6] ] ; 
s [7] : =w [ReadByte (m)] ; s [v [7] ] : =v [s [7] ] ; s [9] : =w [ReadByte(m)] ; s [v [9] ] : =v [s [9] ] ; 
s [10] : =w [ReadByte (m)] ; s [v [10] ] : =v [s [10] ] ; s [12] : =w [ReadByte (m)] ; s [v [12] ] : =v [s [12] ] 
s [14] : =w [ReadByte (m)] ; s [v [14] ] : =v [s [14] ] ; s [16] : =w [ReadByte (m)] ; s [v [16] ] : =v [s [16] ] 
return PermList(s);end;; 


This code will be saved to file chartoperm.txt 
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The function MultEl(el) 


We write elements in generator form as a vector. For example [1,3,2] means 
the element UFL = s[l] * s[2] * s[3] and s was defined above. When we want 
a conversion from the form [1,..,.8] to a permutation with cycles we use the 
function MultEl. Its code can be seen below: 


MultEl:=function(el) 
local x,an; 
an: = () ; 

for X in [1 ..Size(el)] do 

an:=an*s[el [x]]; 

od; 

return an; 
end; ; 


The function Invg(el) 

This function takes an element in generator form as defined above and maps it 
to its inverse also given in generator form. 


lnvg:=function(g) 
local set,x,y; 
set: = [] ;y: = [] ; 

y[l]:=7;y[2]:=8;y[3]:=9;y[4]:=10;y[5]:=ll;y[6]:=12; 
y[7]:=l;y[8]:=2;y[9]:=3;y[10]:=4;y[ll]:=5;y[12]:=6; 
for X in [l..Size(g)] do 
set[Size(g)+l-x]:=y[g[x]]; 
od; 

return(set);end; ; 
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The main program 


The algorithm for the program was explained in the beginning and here we give 
the implementation as GAP code: 


sh:=0;cnt2:=1;positions:=[]; 
for X in sett do 

setl:=MakeSet(10,x); 
set2:=MakeSet(12,x"-l); 
cnt:=0; 

for y in [1..Size(setl [1])] do 

for z in [1..Size(set2 [1])] do 

g:=setl [1] [y] *set2 [1] [z] ; 
if rcv[PermToNr(g)] then 

positions[cnt2]:=Concatenation(setl[2][y],set2[2][z]);cnt2:=cnt2+l; 
set3:=EqClassRel2(g); 

WriteRcv(set3); 
cnt:=cnt+Size(set3); 

fi; 

g:=setl [1] [y]~-l*set2[l][z]"-1; 
if rcv[PermToNr(g)] then 
positions[cnt2]:= 

Concatenation(Invg(setl[2][y]),Invg(set2[2] [z]));cnt2:=cnt2+l; 
set3:=EqClassRel2(g); 

WriteRcv(set3); 
cnt:=cnt+Size(set3); 

fi; 

od; 

od; 

sh:=sh+cnt; 

Print("Positions checked:",sh," ","Positions left:",88179840/2-sh,"\n"); 

od; 

This code will be saved to file mainprg2.txt 
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Running the program and output 


To run the program we make use of all above functions and sets. To read those 
in memory and start the program the following code is used: 


ReadC'makesetv.txt"); #This file contains function MakeSetv 

a:=MakeSetv(l); 

vevenr:=a[1]; 

voddr:=a[2]; 

a:=MakeSetv(2); 

veven:=a[l]; 

vodd:=a[2]; 

Read("makesetvc.txt"); 
a:=MakeSetvc(); 
vevenc:=a[1]; 
voddc:=a [2]; 

Read("makesety.txt"); #This file contains function MakeSety 
yl:=MakeSety(l); 
y2:=MakeSety(2); 

Read("makeset.txt");#Th.is file contains function MakeSet 
ReadC'redf.txt");#This file contains function Redf 
rev:=BlistList([1..88179840],[1..88179840] ); 

Read("permtonr.txt");# This file contains function PermToNr 
Read("eqclassrel2.txt");#This file contains function EqClassRel2 
Read("writerev.txt");#This file contains function WriteRcv 
sett:=Redf(veven); 

ReadC'datav.txt");#This file contains the data included in the paper 

ReadC'chartoperm.txt");#This file contains function CharToPerm 

for X in datav do 

AddSet(sett,CharToPerm(x)); 

od; 

Read("inverseingens.txt");#This file contains function Invg 
Read("mainprg2.txt");#This file contains the Main program 
SaveWorkspace("test2"); 
quit; 

This code is put in a file called start.txt. And all the files(including datav.txt 
at the end of the paper) in the \gap4r4\bin catalog. The program is run under 
Unix with command: 


gap.sh -o lOOOM start.txt 


The program will output: 
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Positions checked:3079007 Positions left:41010913 
Positions checked:18604155 Positions left:25485765 
Positions checked:19391883 Positions left:24698037 
Positions checked:22124037 Positions left:21965883 


And a couple of hours later... 


Positions checked:44089776 
Positions checked:44089824 
Positions checked:44089872 
Positions checked:44089920 
gap> 


Positions 

Positions 

Positions 

Positions 


left:144 
left:96 
left:48 
left: 0 


What does this program do? 

As told in the beginning we need to compute each element in the group Gcube H 
Rc = Ef in at most 22 moves. This program does this. When the program is 
finished the set ’’positions” contains generator expressions of each element in 
Ef. An element of’’positions” has the form [1, 2, 8, 7]. Here 1 means U, 2 means 
E more general n means s[n]. Actually not all elements of Ef are contained in 
positions but only representatives from each equivalence class under relation w. 
To make a final check that confirms that this set contains a representative from 
each equivalence class under relation w the following program can be run: 


Read("multel.txt"); 

corners:=Group(uc,lc,fc,bc,rc,dc); 

rev:=BlistList([1..Size(corners)],[1..Size(corners)]); 

#rcv is a binary vector with all its entries set to true and size 88179840 

for X in positions do 

set:=EqClassRel2(MultEl(x)); 

WriteRcv(set); 
od; 

ent:=0; 

#the program below counts the number of entries in rev that are true 

for X in rev do 

if X then ent:=cnt+l;fi; 

od; 

Print(cnt,"\n"); 

This program takes each element in ’’positions” converts it to a permutation 
with disjoint cycles (this is the form which GAP can handle best). Then build 
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all elements in the equivalence class under relation « and checks them of in 
the bit vector. After this procedure is done the program counts the number 
of entries in the bit vector that has been set to ’’false” and prints them out. 
If it returns 44089920 then every element in A/ has been checked of and we 
are finished. Note that this final check actually proves the whole content of 
this paper is correct. Actually this verification has already been done by Bruce 
Norskog soon after i published the results of this work on the internet. 


Final comments 


It might seem that we were lucky that all elements va Ef had length 22 or less 
but the truth is that we had a backup plan if this would be the case. We find 
generator expressions for each element g in the Rubik cube group in two steps. 
First g is multiplied by 6 G Gcube from the left and then bg G Ef and afterwards 
bg is multiplied by a which takes it to the identity. 

abg = id ^ g = b~^a~^ => g~^ = ab 

Assume that a € Ef has length 24. Assume that the diameter of the Rubik 
cube group is 24. Compute the set 

Xi = {aU~^,aL~^,aF~^,aB~^,aR~^,aD~^,aU,aL,aF,aB,aR,aD}. Then 
each element of this set has length 23 or less. Then a = exi where ei G 
Xi, L(ei) < 23 and Xi is a generator. Then b = X 2 C where X 2 is some generator 
and Cl is an element of length 17 or less. When writing g~^ = ab = eiXiX 2 Ci 
we can choose xi = x^^ and cancellation occurs giving that g~^ has length 40 
or less. It is well known that there exists a position of length 26 so even this 
procedure might fail. So assume that all has length 25 in the list above. Then 
compute the set 

A 2 = {aUU-^,aUL-^,aUF-^, aUB-^,aUR-^, aUD-^,aUU, aUL, 
aUF, aUB, allR, allD} and assume that all elements of this set have length 24. 
Then g~^ = ab = eiXiX 2 Ci. li X 2 ^ U we solve this problem as before. If 
X 2 = U then g~^ = ab = e 2 X\U~^ Ux 3 C 2 . And C 2 has length 16 or less. Again 

we can choose xi = x^^ and the length of g~^ is less than or equal to 40. Fi¬ 
nally some elements in the last set might have length 26 in which case we must 
carry one last step of this. Obviously if there were too many elements in Ef 
having length 24 this computation could take years. I must remark though that 
the elements of length 24 are very rare and very few are known. In 1995-1998 
the only known element of length 24 was the ’’superflip” and today about 150 
others are known to the best of my knowledge. 
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Conclusions 


In this paper we have constructed a program that returns generator expressions 
of length 22 or less for each element in the subgroup of the cube group that 
fixes the edges. As discussed in the abstract this means we have given a upper 
bound of 40 for the cube group. 
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Data vector 

This should be placed in a file called datav.txt: 

datav: = [ ’’aceglpmwqius” , ” acegpkwsujqm” , ” acejhtmoqluw” , ” acejpkwmuqsg” , ’’acekrihownvs” , ” acekvihqmsow” , 

” acelismhqvpw” , ” acelismwohuq” , ” acelismwqhou” , ” aceliwqgonvt” , ” aceliwqugtmo” , ” acelqwmspgui” , ” acepikm- 
sqhuw” , ’’acerikmtohwv” , ’’acesikwmqhpu” , ” acexigqkosmv” , ’’acfgkojqmvxs” , ” acfkihmoqvtw” , ’’acfunilqshow” , 

” acgtikquwnoe” , ” achfikmoqwvt” , ” achliewuonqt” , ” achmikqxosfv” , ” achpiemsqkuw” , ’’achpikfqmsvw” , ’’achrvkmo- 
jsfx”, ” achseuqowkni” , ’’achsikqownfv” , ’’achtikuqonwP , ’’achuislmfxoq” , ’’achvikoqfswn” , ’’achvismoqekw” , ” achviswlonqe” , 
” achxikmvprfs” , ” achxikufmsor” , ’’acjfrnlxosvh” , ” acjgnrloewsu” , ” acjgnrlwesou” , ’’acjgrnewlsou” , ” acjtlnhqevow” , 

” acktieqowhvn” , ’’ackxiegqmsov” , ” ackxiequosgn” , ’’aclgniewprus” , ” aclgniuoexrs” , ” aclgwieqmsup” , ’’aclheipqmvws” , 

” acloniesqgwu” , ” aclonieswguq” , ” aclonisewgqu” , ” aclsvhmqjeow” , ” aclxrimohsfv” , ” acmgeujwoktq” , ” acmgkfjw- 
prus” , ’’acmgkfxqjsup” , ’’acmgpkjwequs” , ’’acmgsljwqvoP , ’’acmoglqfwvjs” , ”acmokfjsqgwu” , ’’acmokfjwqvtg” , ”acm- 
teojwqkug” , ” acmxvgjoqetk” , ’’acnfikhsqpwu” , ’’acnfikhwqspu” , ’’acnfikpxqsvh” , ” acnfikqhwspu” , ” acnfikqxovth” , 

” acnkihvoqetw” , ’’acolihmsweuq” , ” acoqihmewsul” , ’’acoxikmvqsge” , ’’acpiietqmgux” , ” acplismwqveh” , ’’acptieuqmkwg” , 

” acpvikmqesxh” , ’’acqgfnjwlsou” , ” acqokfmsjgwu” , ” acqsikuwmgoe” , ’’acruikomxsge” , ” acsuikomgeqw” , ’’acsvri- 
hownke” , ” acsxikfoqngv” , ’’actfikmqvhow” , ’’actfikmuqgxp” , ” actfkmjwqhou” , ” actlvieqmhow” , ’’actpieuqmkgw” , 

”acugfnjwoslq” , ”acugkfojwsmq” , ”acuokfmswgjq” , ”acuokfsmjgwq” , ’’acvfikmswphq” , ’’acvfqjmkpsxh” , ’’acvgpetk- 
mjwq” , ’’acvlhjmoqftw” , ”acvmietoqkxh” , ”acvmikqfosxh” , ’’acvnkfxqjtgo” , ”acvriktomewh” , "acvsilmoqexh” , ’’acvwikm- 
rosge” , ’’acwtieqgokmv” , ’’acwuiemoqslh” , ’’acwuiksqmego” , ” acxfikhmupqs” , ” acxfikoqmsvh” , ’’acxlisumfhoq” , ’’adegmki- 
wpqus”, ” adegqkmupitx” , ” adeohkmwqius” , ” adesukhomiwq” , ’’adeuokqsximg” , ” adeuqkmgtiow” , ’’adewqomglius” , 

” adexikmsqoug” , ” adfvslnohiwq” , ” adgfqkmvxitp” , ” adhuisrwpnle” , ” adkgmsfwpivq” , ” adkgqmewpius” , ” adkgws- 
fqmipv” , ” adknqgpuwifs” , ” adksowhfuiqn” , ” adktiomrhvfw” , ” admfkuhxrisp” , ” adpgikxeunqs” , ” adpgslmwrifv” , 
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” adpgslwfuirn” , ” adptkgnwqifv” , ” adrxsnlohifv” , ” adsvowqlhime” , ” advlmerwpgti” , ” advsiemkqoxh” , ” adxkmeuqpgti” , 


” aegoidmlqsuw” , ” aehtikmcquwo” , ” aehtikmvqwdo” , ” aekgismvqwdo” , ” aelgniwcuqso” , ” aesoiwmvqclg” , ” aesuqk- 
mgpwdi” , ” aevdikmhqwto” , ” aexoivtgmqlc” , ” afckiumoqgtw” , ” afkoismvqwgc” , ” afkrcsmxipug” , ” aflcqhpuwims” , 

” aflcqhupmiws” , ” afqckhjuwpms” , ” agclnpewqjus” , ’’agcuikmoqswe” , ” agdsiemoqluw” , ’’agdsikfqmxov” , ’’agecltqow- 
imu” , ” agecqkmpxius” , ” agecskmoqiuw” , ’’ageisdmoqluw” , ’’agejckmwqpus” , ’’agejktmoqcuw” , ” agekciurmpws” , ”agel- 
cjqosumw” , ” agelcjurmspw” , ” agelcpmrixus” , ” agelipquwcms” , ” ageliqmwocut” , ’’agelismoqvdx” , ” agenckuoitwq” , 

” agenckuqitow” , ” agenikqxocvt” , ” agenikuwqcot” , ” agenksrdpiuw” , ” ageoikmdqtuw” , ” ageoikmrxsuc” , ” ageoikm- 
sqwuc” , ”ageoikmsqxud” , ’’ageoikmsxcur” , ’’ageoitmlqcuw” , ’’ageonkqcjsuw” , ”agepltmricuw” , ”ageqckuniswp” , ”age- 
qxkmcosuj”, ” agesckuomxri” , ” agesikmcqwuo” , ’’agesikmrxcuo” , ’’agesldmoqiuw” , "agesnkuwocjq” , ” agesqkmwcoui” , 

”agesqkpuwidm” , ”agesqkuomidx” , ”agesrkjomcuw” , ”agesrkjwmcou” , ”agesvkjwocmq” , ’’agesvkojmcqw” , ’’agesxk- 
moucqj” , ” ageuckiosxmr” , ’’ageuckqoixmt” , ” ageuikmwqsoc” , ” agevismoqlxd” , ” agewikmsoquc” , "agewikmsqvpc” , 

” agewikmvqsdo” , ” agewikqucsmo” , ’’agewikucmsoq” , ’’agewikuqmtdo” , ” agewmkuordsi” , ” agewpkmjcsuq” , ’’agexikm- 
cpsuq” , ”agexikmcqpus” , ”agfkvqxdmitp” , ”agflniquoctw” , ”agfsnivkqcpw” , ”agjsfvmwoclq” , "agknwsjoqcfv” , ’’agkrc- 
smoivfx” , ” agksiqfomcwv” , ” agksiwmqecvp” , ” agksrwmfocjv” , ” agktipmwqcfv” , ” agkuiwtoqcne” , ” agkvpsfncqwi” , 

”agkvqsmoxifd” , ’’agkvsfjoqcnw” , ”agkvwsmpqidP , ”agkwiomcqsfv” , "agkwiomrtcfv” , "agkwiomsqfdu” , ’’agkwismc- 
qvfo” , ’’agkxcjmoqsfv” , ’’aglcrtmoweju” , ” aglfixuomctq” , ” aglovieqmscw” , ’’agloximcquse” , ” aglrniewpcus” , ’’aglsjweo- 
qcmu” , ’’aglwvieomqsc” , ” agmctkjoqeuw” , ’’agmldpjwqeus” , ’’agmodejvqktw” , ” agmsikqoecuw” , ’’agmsikqwecou” , 

” agmsikuwoceq” , ’’agmwkfjcqsuo” , ’’agnfiwquoctl” , ’’agnrsljowceu” , ” agodikmvqswe” , ” agplismwqcfv” , ’’agpsikfwqcmu” , 

”agqsikeomcuw” , ”agqsikewmcou” , ”agscikmoqeuw” , ’’agsjrkeowcmu” , ”agswnielqcuo” , ”agtdikmwprfv” , ’’agtovqed- 
mikx” , "aguerisowcml” , "aguliemoqcxs” , ’’agusikewocmq” , ” agusikfdpnwq” , "aguwkfjqmsoc” , ” agvcelmoqjtw” , ’’agvcflmo- 
qitw” , ’’agvcfpmiqktw” , ’’agvfikmcqspx” , ’’agvfkmjoqctw” , ” agvkiemcqsxo” , ” agvkiemwqods” , ’’agvrcemoiltw” , ’’agvrce- 
moitxk” , ” agvsipmxqcle” , ” agvwiemcqkto” , ” agwsipmlqcue” , ” agwvieqtocmk” , ” agxkcetnupqi” , ” ahecitmrkpuw” , 

” aheodirsxnul” , ” ahesqkmowdui” , ” ahnuieqcowlt” , ’’aidkmersowhv” , ” aiegckqwnpus” , ” aiegquptwcmk” , ’’aiehmutq- 
pldx” , ” aiekqumdpsxh” , ’’aielqspdwnvh” , ’’aielqspudhmw” , ’’aienckwguqsp” , "aienpkwcuqsg” , ’’aieomkugqwsc” , ”aier- 
ckmsgwup” , ” aierpkmscwug” , ’’aieskhmdpruw” , ’’aiesmkuwqgoc” , ” aiespkmdgquw” , ” aievqkpdwhms” , ” aiexckmgqusp” , 

” aikgcsmqvpfx” , ” aikhscqoxnfu” , ” ailgncesqxup” , ” aimgekupqsdw” , ” aimgslpfdvrw” , ” aiswceukmroh” , ” aiurewd- 
kmptg” , ’’aiurkfdgmwsp” , ’’aivgcewktnqp” , ’’aivrcehkmwsp” , ” aixfmkptdhuq” , ’’ajchkfosxqnv” , ” ajegpkuwcmsq” , 

” ajegxkqtcpmu” , ” ajeovkmwcqsg” , ” ajesqkpdhmuw” , ’’ajewtkucmqgo” , ” ajkgxsmrcpfv” , ’’ajkvmsfdhqwp” , ’’ajndwep- 
kvqtg” , ’’ajogmewcuqtk” , ” ajsgpvqlcwme” , ” akecqhiwmpus” , ’’akecwhirosmu” , ” akehctmoiquw” , ’’akeoiguqmwsc” , 

” akeoiqncwshu” , ” akerchmsixup” , ’’akeritmwphuc” , ” akeschmriouw” , ” akesiwugmcoq” , ” akhpcsmwiqfv” , ’’akhtqmew- 
pjdv” , ” akmchfjrosuw” , ” akmuetjcqwph” , ” akqhcemsiwup” , ” akvcmetiqpwh” , ” akvncetgiqwp” , ” akvscemriphw” , 

” akvsiehqmcow” , ” akvtweqgpimd” , ” akvwceitgpmr” , ’’akxoivtemqhc” , ” akxqcemgipus” , ’’alcgifmwqous” , ’’alcoifug- 
mqsw” , "alcsrfmowgju” , ”alcsvfmojgwq” , ”alcsvfmqjgow” , ”alegctirwpmu” , ”alegmtudpiwq” , ’’aleqptmjgcuw” , ”ales- 
diqowgmu” , ” alesdiuomgwq” , ’’alesixmoqcug” , ’’alesoidrnguw” , ” alhsiemoqcuw” , ” almehtrdpiuw” , ” almgqsjwpdfv” , 

” almsfdjoqguw” , ’’alpscfmriguw” , ” alsgivmwqoec” , ” alsgrieowcmu” , ’’alsgvieomcwq” , ” alvntfjoqcxh” , ” alvsdemo- 
qgxi” , ” amcgikouwqes” , ” amegiksqcuwo” , ” amegikuqcwso” , ” ameoikwruchs” , ” amepqkwguisd” , ” amesikgocquw” , 

” amesikocwuqg” , ” amesikoqugwc” , "ameuhkscjqwo” , ” amlgipewqcus” , ” amrfikuhscwp” , ” amswpvceqjlg” , ’’amugk- 
fwcjqso”, ” andtgvlqipfx” , ” anegikqcwups” , ’’anigkfquwpds” , ’’anksiwqhdufo” , ’’ankwisfgqcuo” , ’’anlgriewpcus” , ” anlpchuqiews” , 
” anlvgqepxsdi” , ’’anpuiglcqsfw” , ” anvkjcewqpth” , ’’anvwcepqitgk” , ’’aoclrfwjugms” , ” aocuikqswgme” , ” aoecikquwtmh” , 

”aoecikuqmtwh” , ”aoecktmwqguj” , ’’aoednskwqjuh” , ”aoedqkmvhitw” , ’’aoedukwirgsm” , ”aoegdmjwqkus” , ’’aoegikm- 
rdvxs” , ” aoegikqsucmw” , ” aoegikqtvcmw” , ” aoegikquwntc” , ” aoegiktqucwn” , ” aoegikuxmctq” , ” aoegikwqmctv” , 

” aoegisukmcwq” , ’’aoegltmwqiuc” , ” aoegnkjqucws” , ” aoegskmqjcux” , ’’aoehskmdqiuw” , ” aoejmkdwqgus” , ” aoejrk- 
suwcgm” , ’’aoelisdqmhwv” , ” aoelismdqhxu” , ” aoelituqmcwg” , ’’aoelnrcuwgjs” , ” aoelqwmvhitc” , ’’aoenkgj vqctw” , 

” aoenskwjvcqh” , ” aoerckmwigut” , ” aoerckmwjgus” , ’’aoerckwmuhis” , ’’aoerclmwigus” , ” aoerjkdwngus” , ’’aoerskmx- 
icug” , ” aoesihmwqkud” , ” aoesihquwcml” , ” aoesikmcquwg” , ” aoesukwrhcim” , ” aoetikmgqcwv” , ” aoetikqhwcmu” , 

” aoevhkjwqcns” , ” aofgimuwqctk” , ’’aofgksquwcni” , ” aoflsdjwqgnv” , ’’aofriktwmcuh” , ” aogcrnlwfius” , ” aogmvilqecws” , 
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” aogrldewnius” , ” aohckfjwqnus” , ’’aohclsqfwimv” , ” aohlivmwqcfs” , ’’aohlnevwqjds” , ’’aohlvemqjcws” , ” aohlweiuqcms” , 

” aohtifmwqcul” , ” aohtikfdmvrw” , ’’aohwikmuqcfs” , ’’aoignsquwckf” , ” aoigrkewncus” , ” aoilnevwqchs” , ” aojgfnlsxruc” , 

” aojgkvcqmews” , ’’aojglrwfuctn” , ” aojgnveqlcws” , ” aojkngewqcus” , ” aojtndkwqgfv” , ” aojwfnqcguls” , ’’aojwgvsqm- 
cke” , ’’aokdqshfwimv” , ” aokgcsirvxmP , ’’aokgismvwcfq” , ” aokgiswufcmr” , ” aokgixmuqcfs” , ” aokgmsqwicfv” , ”aok- 
grsjwmcfv”, ” aokgvsjfqcmx” , ” aokgxsqfucjn” , ’’aokjdemwqgus” , ’’aoknigfwqcus” , ’’aokrsgjfwcmv” , ” aoksivmwqcfg” , 

”aokuismhqcfw” , ’’aokvishfrcmw” , ”aokwiemgqctv” , ’’aokwismgqufd” , ’’aolcniewqtuh” , ”aolgnvjwqces” , ’’aolgwieqn- 
cus” , ” aolhismwqcfv” , ’’aolmeiswqcgv” , ’’aolsnihurcew” , ” aoltiwhurcfn” , ’’aoluiqswncge” , ’’aomchvlqejws” , ” aom- 
ckfhwqius” , ’’aomckfjwqtuh” , ” aomedkgwqjus” , ” aomedkjwqtuh” , ’’aomgikeuwcqs” , ” aomgxtjlqcue” , ’’aomidfhwqlus” , 

” aomligrwecus” , ’’aomtfljwqgdv” , ” aomtilcwqgfv” , ” aonfikdsqhwu” , ” aonfikqxhcvt” , ’’aonfwkqtucih” , "aonlisqfwgdu” , 

” aonrceiuwhtk” , ” aonuiewckgqt” , ” aonuisfwqcgk” , ’’aonxieqkvctg” , ’’aoqcehjuwkms” , ” aoqgntwlucje” , ” aoqhniesx- 
ckv” , ”aoqhskjfwcmv” , ’’aoquismlwche” , ”aosgrvmlwcje” , ’’aosgvxmqjcle” , ”aotfikmdwhuq” , ’’aotliemdqhuw” , ”ao- 
tungkwqcie” , ” aougfrljmcws” , ” aougriwlfcsn” , ” aoulcfnqjgws” , ” aounkfdqihws” , ’’aoutigqlwcme” , ” aovcgemwqitk” , 

” aovciemwqshk” , ” aovcikmhqsxf” , ” aovfikmhqcwt” , ” aovfikmwqtgd” , ’’aovgemjdqktw” , ” aovliemtqcxg” , ” aovrceitwhmk” , 

” aovsikmdqgxf” , ’’aovtikmgqcfx” , ” aowckhmeujqs” , ’’aowgnielucqs” , ” aowjkfmuqgsc” , ” aowliemsqcgv” , ” aowuieqt- 
glmc” , ”aoxcikfqmtuh” , "aoxgniesqcuk” , ”aoxrceuimhtk” , ”apegikqwsnuc” , ’’apegmkuqxisc” , ”apejmkugxqsd” , ’’apels- 
dmriguw” , ” apevckqxihms” , ” apknisfvqcwg” , ” apvnkfdgiwqt” , ” aqcwikmfsogv” , ” aqegcomixlus” , ” aqegikmdsowv” , 

” aqegvkmcjwso” , ” aqelismdxouh” , ” aqencksuwpig” , ” aqeoikmgcsuw” , ” aqeoikutmwdg” , ” aqesvkmwjgoc” , ’’aqevikxs- 
modh” , ’’aqewikucmosg” , ”aqewvkmgjcso” , ’’aqjsfnuwlgoc” , ’’aqkgosmiwufd” , ’’aqkncgutixfo” , ”aqlwniuoesgc” , ”aqnf- 
pkcjwsvh” , ’’aqtvkfjgmwdp” , ” arcgienuwptk” , ” aregxkumdpis” , ” areldinuwphs” , ’’arflgsmwipdv” , ’’arhojckusxfn” , 

” arlgcinuwpes” , ” artocvhimxle” , ” arvwceipgtmk” , ’’asedikupmqwg” , ’’asegikwpudnr” , ” asegipulmqwc” , ’’asegixqul- 

cmo” , ’’asegvkmqjwco” , ” asekigmoqcuw” , ” asekiumgqwdp” , ’’aselcwmgipqu” , ” aselihmoqcuw” , ’’asenckigqxuo” , ” aseoikuqgmwc” , 

”aseoikuxmcgr” , "aseqikmogcuw” , ”aseqipulmcwg” , ”aseqpkujhmwd” , ’’aserckmviwog” , "aseuikqcwomg” , ’’asevikmo- 

qchw” , ’’asewcqmlipug” , ” asewikmogquc” , ’’ashxikfoqncu” , ’’askoigfwqunc” , ’’aslcngeiqwuo” , ’’aslorimweugc” , ’’aslqvimwgceo” , 

” aslrniegxouc” , ” aslwipmgqeuc” , ’’aslwnieohcur” , ” asnuiepwqcgk” , ’’asqofnlugwjc” , ’’asuwikfqmcpg” , ” asvkiemgqwpc” , 

” asvwiemkqodg” , ” aswvioqlgcme” , ” asxhikmcqofv” , ” ategikmvqowc” , ” atelcuqhiomx” , ” ateliwmchruo” , ” ateqhk- 
mociuw” , ” aterxkigcvmo” , ’’atngceikuxqp” , ’’atnuweqkpdjh” , ’’atpgikudmqfw” , ” atvgpeikcxmq” , ” atvwiemkqodh” , 

” auegckqswnip” , ” auegikmxqotc” , ” auegiknswrdp” , ’’auejpkhwdqms” , ’’aueoikqgtxmc” , ” aueoqkmsxdgi” , ” aueoskigqwmc” , 

” auesikmwqocg” , ’’auesrkiwhomc” , ’’aufmriskwodh” , ’’auhtwkinqofc” , ” aulsniqogcew” , ” avcfikmsqwph” , ’’avegckm- 
siwor” , ” avegikmcqspw” , ” avegikmsqodw” , ” avegikomqsdw” , ” avegikqdswmo” , ” aveligmoqsdx” , ” aveoikmcqsxg” , 

” avesikmwqodg” , ” avesqkmdpwjg” , ” avkocgfmixqt” , ” avkpismhqwfc” , ’’avkwismohcfr” , ” avkwmsfdhioq” , ’’avwsiemgqcpk” , 

”awdjlhmsqofv” , ”awedismkqouh” , ’’awedqgpukims” , ”awedqkphvims” , ”awegikmdqosv” , ’’awegikmuqodt” , ’’awegikom- 
suqc” , ” awegikqscumo” , ” awegmkurpcsi” , ’’awegqkmipvtd” , ’’awegskqupdmi” , ” awejckosnqug” , ’’awekisuqmodh” , 

” awelivmhqods” , ” awelpdmsirug” , ” awenckvsiroh” , ” aweoihmsqluc” , ” awerckusmioh” , ” awesikuomgqc” , ” aweuik- 
moqsgc” , ’’aweuikmsqcgo” , ” awgjpvmscqle” , ’’awjgfvmsqolc” , ” awjrcesmuphk” , ’’awkoqsmduifg” , ’’awkoqsmvhifd” , 

” awmfikhvqods” , ” awmoekjdqgus” , ” awmokfjgqusc” , ” awmuktrepdgi” , ” awpfikmuqcth” , ” awrfikdusomh” , ” awse- 
hvmlqcjo” , ’’awsmikgeqouc” , ’’awsoeqjmuchk” , ” awvoremtkcjg” , ’’awvsikfqmodh” , ” axegikusmqdo” , ’’axeoqkmsg- 
dui” , ” axhtikqfucpn” , ” axksiqocfumg” , "axphckifsvqn” , ’’axqgemitupdk” , ” bdmuqtlxjoeg” , ’’bfmsokjhqcuw” , ’’bkeogjmd- 
sxuq” , ’’bwehvqkpjdnt” , ’’caljnpegqsuw” , ’’camlxfroigus” , ” cbsgivlqmpwe” , ” cbvfqjmwplgs” , ” cealnigousqw” , ’’ceglib- 
mvqwps” , ’’cejlnpawqgus” , ” cesagjuwmqlo” , ’’ceslibmguwqo” , ’’cfsulranwoig” , ” cfsuqjmbphxk” , ” cgalfiuoswmq” , 

” cgaliqewosun” , ” cgalpiuemsqw” , ” cgalqiemuswp” , ” cgaltiqouemw” , ” cgaltiuomeqw” , ” cgameiwoqslv” , ” cgateiuq- 
know” , ’’cgauniswoelq” , ’’cgaurimwseol” , ” cgbxkimoqsfv” , ’’cgeiqpmbxlus” , ” cgelibmvprtw” , ” cgelibqtowmv” , ’’cgelib- 
vqmsox” , ” cgelibwusqom” , ” cgelmjpwbsuq” , ’’cgelniuwosaq” , ” cgelpbujmsqw” , ” cgelqbouismw” , ’’cgelwbiornus” , 

” cgesqamokiuw” , ’’cgetibmlovwq” , ’’cgfbniqosuxk” , ’’cgfxblniqsov” , ” cgfxikmoqsbv” , ’’cgjlbfqwmsou” , ” cgjlfsaoqxnv” , 

” cgkasjqowvfn” , ” cgktibvmfqow” , ’’cgmakpjwqeus” , ” cgmesbjoqlxu” , ” cgmilfrotauw” , ’’cgmpbljrtwfv” , ’’cgmxufao- 
qkis” , ” cgnlibfowrvs” , ” cgnlibfqusow” , ” cgofibqmxsvk” , ” cgpilnfqavxs” , ” cgplibufmsqw” , ” cgplibuqfsmw” , ” cgp- 
saeuqmixk” , ’’cgpxikfqmsub” , ” cgrfibwuokmt” , ’’cgsfibmuqkow” , ’’cgsubfjwomlq” , ” cgsuibqfomkw” , ” cgtbniuoeqxk” , 
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” cgulibfnqsow” , ” cgvbiemoqlxs” , ” cgvlbfrojnxs” , ’’cgvnafrokixt” , ” cgvsiembqokw” , ” cgvtibwoqkfm” , ” cgvxniqboetl” , 

” cgwlbfqjosmu” , ’’cgwlibeqosum” , ’’cgwlibqeusom” , ’’cgwlifmoqvta” , ’’cgwtblqfjvom” , ’’cgxfnbutokjq” , "cgxtibmo- 
qlfv” , ” chelibpuwnrs” , ” chesikunrwob” , ’’chvnqfptwbjk” , ” chxlqemtpiub” , ’’cimalfrwpgus” , ’’cixkemugpqtb” , ”ck- 
afnisgquwo” , ’’ckafsiomvwgq” , ” ckahniuseqwo” , ’’ckgbniseqwpv” , ” ckwhibmoqvte” , ” ckxpieuqmath” , ’’claovimgsweq” , 

” clebiomgqsuw” , ” clergbmoisux” , ’’clewgbmrisup” , ” clewibmogsuq” , ” clmgbfusjwoq” , ” clmsbfjwqgou” , ” clqajfhu- 
osmw” , ’’clqsbfmwjgou” , ’’cluwbfjqmsgo” , ’’cmalriseguwo” , ’’cmulbfjqgwso” , ’’cnebiqkpvgws” , ’’coafnlrwigus” , ”coaf- 
sjmwqguk”, ” coaheikmuxqs” , ” coalnwiuqges” , ’’coalriesugmw” , ” cobgilqfwsmv” , ’’coegikuamqws” , ’’coegmkuaqiws” , 

” coehibmuqxtk” , ” coehibukrnws” , ’’coehqbkuwims” , ” coeiahmwlqus” , ’’coelrbmsvxjg” , ” coelsbujmgqx” , ” coergaiuwlms” , 

” coernbtjvgwk” , ” coesibmvqlhw” , ’’coetbgjwqlnv” , ” coetkbxqjgnv” , ” coglbtjqmuwe” , ” cohtbljqmuex” , ” cohuibfqmwtk” , 

” coifrbvwngtk” , ” cojhenaqkvws” , ” cokfibnqsgwv” , ’’cokithewnrbv” , ’’cokjnaquwges” , ” cokrvbmxjgfs” , ’’coktibfhuxnq” , 

” coktibfqmxug” , ” coljnberhsuw” , ” comhbljfwtuq” , ’’comlfubqjgws” , ’’conbrixutgle” , ’’confaktwqiuh” , ’’confibtkquwg” , 

” confibuqwksh” , ’’conwibqfugtk” , ’’coqlbfwjmgsu” , ” coqlfnsjwgau” , ’’coqlibesmgwu” , ’’coqlnbvwjgfs” , ’’coqlviweagms” , 

” coqtbfjukgmw” , ” coqwbtjufgln” , ” corlibwmfvtg” , ’’corlibwusneg” , ’’cortliawfgnv” , ’’cosfibqmwguk” , ” cosnbfjuqgwk” , 

” cosqvnweaglj” , ’’cosubkmwqfgi” , ” coulbfsjqgmw” , ” coulibwmqges” , ’’covfibhwnktq” , ” covfmbhwrkti” , ” covfnbkwqgjt” , 

” covgnaewqski” , "covlrbjwegns” , ’’covwisfqmglb” , ” coxfibumtgqk” , "coxknbteugqi” , ’’coxtiulqmgbe” , ” cqaoxiemugsl” , 

” cqaxnikovgfs” , ” cqehmbwpuisl” , ’’cqelibuspgwn” , ” cqewihusmola” , ” cqjbenxskoug” , ” cqmlsbvxjofg” , ” cqsoibml- 
wgue” , ”csagniuoeqwl” , ’’csahnqelpwui” , ’’csalniouwqge” , ”csaqnikohwfv” , ”csawjkegqnuo” , ’’csfbnipgqwuk” , ’’csnhibe- 
owkru” , ”csqlubigwemo” , ’’csvfrbmxgkjo” , ’’csvmikqgwfao” , ”ctbrgoixlnue” , ”ctelawmiqoug” , ’’ctvgweikqbmp” , ’’cuakrip- 
swgfn” , ’’cueagsmpirxl” , ” cuelibmgxqto” , ” cukgibqnosfx”, ” cunfaoiswrhk” , ’’cvbrgsmliwfp” , ” cvksibmoqwfh” , ”cvk- 
tfiqgmwbo” , ’’cvmpbtlqjwge” , ” cwehbujsqkno” , ” cwelmbusjgor” , ’’cwelqbumiosg” , ” cwevnbikqotg” , ” cwjlnveqaosg” , 

” cwkhbsjqmoev” , "cwluiomsqhae” , ’’cwpakimsqhfv” , ” cwqpbljfsvmg” , ” cwvhifqbkoms” , ’’cwvnbfjrpgtk” , ’’cwvnofrsjlhb” , 

”dbegikuqmows” , "ebmgksrwpiud” , ”ecmlsajoqguw” , ”ecpnbkwjusrh” , ’’ecpnbkxjusqh” , ”ecsxbknqjvgp” , ’’edmkqb- 
hwpius” , ” edphknasxiqv” , ” egloibmcqwsu” , ’’egsaidwmouql” , ” egsvcojrbnlw” , ” egvkbmitrpdw” , ’’ehpkqbjsxdnv” , 

” ejulmbrpgdws” , ” ejusbkgocmwq” , ’’ejusbkndhqwp” , ” eknscqjhbovx” , "ekvgbmcqwtoj” , ’’elaoidugmqsw” , ’’enkobgqvswjc” , 
” eogcxkumjaqs” , ” eokgcsjrbxnv” , "eomgjkdwqaus” , ” epksbwduihrn” , ’’epvgbmtkrdwi” , ’’erhucsblwpjm” , ’’erxgb- 
vtidplm” , ” esuwbkmogcjq” , ” euigbkxsqomc” , ’’euinbkqgswdp” , ” euinbkqtogxc” , ” evajpdwnglqs” , ’’ewctmibvkqoh” , 

” ewmockjhbrus” , ’’fcagsqjolnuw” , ’’fhspbmdviwlr” , ” fhvcqjwmtapk” , ” fonaidqkwshu” , ” gcajntuqelow” , ’’gcalitmo- 
qeuw” , ”gclantuoejwq” , ”gdeaqjmswlup” , ’’gemxdaqvoljt” , ”gfurbdjomksx” , ”gipkecmqbtwv” , ’’gkanqdewpjus” , ”gk- 
canqpuwies” , ’’gkeaspquximc” , ’’glcaxiouesmq” , ” glebmurdpsxi” , ’’glmrdftojbux” , ’’glnfidqbwspu” , ’’glpxvcbimsfq” , 

” gnparceuwlis” , ” gobwidkqmtfu” , ’’gockrbjwmtfv” , ” gojsdvaqmwle” , ” goqjnbxudles” , ” gosurdmwelja” , ’’govlidmwqbtP , 
’’gpeiqdusmbwl” , ’’gquednaljwso” , ” gsajvdolewqm” , ” gseiwcqapuml” , ” gtdkqjmwpafv” , ” gtnuidqeolxb” , ’’guitdbqn- 
wofl” , ”guktidfqmbow” , ’’gwbtmjpvqdfl” , ’’gwecipquslma” , ”gwvfidmlqsao” , ’’gxplrcnfitub” , ”hcofibqmusxk” , ”hra- 
jsdfnkxov” , ” ibfogqvxmcsk” , ’’ibhoqkmvxcft” , ” ibpgqkfdmvxs” , ” ibptkgnwqcfv” , ’’icbhqwmsplfv” , ” icelqwmgosub” , 

” icvomehwrksb” , ’’idaukmsqheow” , ’’idegqkusmowb” , ” idetqlwmbgpv” , ” idkhbfrnxout” , ” idsulrafpgwn” , ” idvoqextbgml” , 

” idvtqblwpgfn” , ” iektqbmvpwhd” , ’’iekwbvgscorn” , ” ifvkqbmsxodh” , ” igqnblpfdsuw” , ’’ihbksdwfvoqn” , ’’ihclnqpuwfbs” , 

” ihpbnqukedxs” , ’’iksuqemhxbdo” , ’’ilbxqdhuofms” , ’’ilvbqemdxohs” , ” imrfwbuhsdpk” , ’’iodglbmwqtfv” , ’’iomfbutxrgdk” , 

” ipelgbusqnwc” , ” irsbdhlvmeow” , ’’isalnqdguexo” , ” isepqbxukdmg” , ” iuehqbwmkdps” , ” ivatnqhfdwko” , ’’ivbpnqewkgdt” , 

” iwtfqbmkpuhd” , ” jnakwuqdsgpe” , ” joqgebvwmctk” , ” jragvkpwmces” , ” jskwnqugecbp” , ” jsowndegqlua” , ’’juaonetwqchk” , 
” kdjebnrohsuw” , ’’kdqgwsanpfjv” , ’’kesugdinoqaw” , ” kgiebdmoqvtw” , ” khmaqdjofsuw” , ” kibtgfnqdvow” , ’’kixfqbupm- 
cth” , ’’krmxcsjobgfv” , ” ksmwcuafhqpj” , ” Iceoshiaqumw” , ’’lecunipaqshw” , ’’licmgovbfwtr” , ’’Inadqtejhwuo” , "Iphvc- 
sjrbenw” , ” Iwepchsmuria” , ’’Iwjbhfrspdnv” , ’’majedvoqglws” , ’’mgarcoelisuw” , ’’mhqlvawtepjc” , ’’mjeskbuqcgow” , 

” mjpfrbukcgxs” , ” mleapdugsqwj” , ” mocwblvqjhte” , ’’mougkiawqces” , ’’mqoclsfjavgw” , ” mruagfjpdlws” , ’’mrugfk- 
wjdpbs”, ’’ngibqdxlsufo” , ’’obxngvtciqkf” , ’’ojdkxrmsgafv” , ’’olfbndquwgti” , ” omriwdtbkvgP , ” pdmlwfargsju” , ”pg- 
bicnqxfsvk” , ” prjncuawfkhs” , ’’psmwcfjlbgur” , ” pvkacimgswfr” , ” qameowucjlsg” , ’’qjecwsmgoaul” , ” qkenixvocsha” , 

” qwvfcbmoikhs” , ’’rnoticvelxag” , ” scegaomrixuk” , ’’sgelcbqowimu” , ” sgwoniuqelac” , ’’tgeawdmoilur” , ’’trgiwpmbe- 
clv” , ” twvoiemhqadk” , ’’txkpbfdgiurn” , ” unkgasfqiwco” , ” utbwiemhqocl” , ”vjtkoqmbhcex” , ”voekbsjwqcng” , ’’wjeovk- 
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maqcsg” , ’’wseclpmgaruj” , ’’xoahcqesniul” , "xsvnafjgkrdp” ]; 


The last data have been selected from a list that would take about 2 months to 
generate on an usual PC but thanks to my advisor Gert Almkvist access to fast 
computers have been given and made it possible to compute in just a couple of 
days. 
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